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ABSTRACT 


The dynamics of a second-order phase transition during inflation, which is induced by time- 
variation of spacetime curvature, is studied as a natural mechanism to produce topological defects 
of typical grand unification scales such as cosmic strings or global textures. It is shown that their 
distribution is almost scale-invariant with small- and large-scale cutoffs. Also discussed is how 
these cutoffs are given. _________ 
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1 Introduction 


< 


The idea of spontaneous symmetry breaking plays an important role in constructing unified the- 
ories of elementary interactions. At ultrahigh temperatures of order of their unification scale, it 
is believed that the presently broken symmetries are restored due to high-temperature correction 
to the potential of Higgs fields which determines the symmetry of the system. Unfortunately, 
however, it is hardly possible to verify their predictions experimentally, since the energy scale is 
higher than the accessible scales by many orders of magnitude in a laboratory. 

Trying to get out of such a frustrating situation, many people applied these theories to the 
early universe whose temperature, according to the conventional big bang cosmology, was once so 
high as the unification energy scale. One of their natural consequences is that the universe has 
presumably undergone a number of thermal phase transitions in the course of its early evolution, 
in some of which topological defects such as walls, strings, monopoles, or textures may have been 
produced through the so-called Kibble mechanism [1]. Among them cosmic strings or global tex- 
tures may help large-scale structure formation if their energy scale is O(10 16 GeV) [2, 3]. Detection 
of magnetic monopoles, on the other hand, is an important clue to verify unified theories because 
they are inevitably produced when the electromagnetic 17(1) gauge group branches off at some 
symmetry breaking [1, 2, 4]. 

In discussing these phase transitions, it has been implicitly assumed that the universe was in 
thermal equilibrium state at least by the GUT era. However, in order to attain thermal equilibrium 
from an arbitrary initial state of the universe, it is necessary that particle interaction rates exceed 
the cosmic expansion rate, which is not always possible in such an early stage of cosmic evolution 
[5]. Hence formation of topological defects with typical GUT scale may not be described by the 
Kibble mechanism correctly unless the universe started its evolution in a thermal equilibrium 
state. 

It is more natural to expect that our universe started its classical evolution out of a chaotic state 
governed by quantum and thermal fluctuations [6] and that it underwent inflation to be globally 
homogeneous and isotropic as observed today [7]. Then it turns out to be after the reheating 
epoch that the universe was first filled with radiation in thermal equilibrium. Unfortunately, the 
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maximum temperature it experienced, or the reheat temperature, may not be so high as the GUT 
scale generally in order to avoid too much gravitational waves or density fluctuations [8] and/or 
too many gravitinos [9] to be produced after inflation. On the other hand, models with high 
enough reheat temperature or extended inflation scenario may not solve the monopole problem, 
even if they may keep large enough density of strings or textures. Thus it is very difficult to obtain 
a sensible scenario of the early universe which provides an appropriate initial condition of galaxy 
formation through topological defects in the grand unification scale. 

In order to resolve this difficulty several mechanisms of non-thermal phase transitions have 
been proposed in which the Higgs field is coupled either with spacetime curvature 72 [10, 11] or 
the inflation- driving field <p [12, 13, 14]. The former mechanism is especially plausible since the 
effective potential naturally has a finite-curvature correction in the inflating spacetime just as it 
would have a finite-temperature correction in the hot big bang universe [15]. In these scenarios 
phase transition takes place during the inflationary stage due to time variation of 72 or <£ and it 
is triggered by quantum fluctuations rather than thermal fluctuations. 

While the properties of thermal phase transitions have been extensively studied, those of the 
above non-thermal phase transitions have not been fully investigated. The purpose of the present 
paper is to clarify the dynamics of a second-order phase transition during inflation as well as 
the spectrum of the defects produced. Complementary to the present paper, the case topological 
defects are produced through a first-order phase transition has been studied by Copeland, Kolb, 
and Liddle [16] in the context of extended inflation scenario and quantum creation of defects 
through tunneling during inflation has been investigated by Basu, Guth, and Vilenkin [17]. 

The rest of the paper is organized as follows. After introducing our model in §2, we examine 
when evolution of the scalar field becomes deterministic in §3. Then we estimate the power 
spectrum of quantum fluctuations in §4. In §5 the result of numerical simulations on the spectrum 
of defects thus produced is reported and its analytic interpretation is given. Finally §6 is devoted 
to discussion and conclusions. 
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2 The model 


In the present paper, for definiteness, we consider the chaotic inflation scenario realized by a 
massive scalar field <f) with mass m in the spatially flat Friedmann-Robertson- Walker spacetime 
da 3 = dt 3 — a(t) 3 dx 3 = a(ji) 3 (dT) 3 — dx a ). In the inflationary stage, the classical evolution of <£(£), 
the scale factor a(t), and the Hubble parameter H(t), is given by 


4>{t) 

= 4 

mM pl ( . . V 

0 2VS£ ( ‘ ‘ o) ’ 

(1) 

m 


47rm 2 , . m 3 , 

(2) 

= ^ 

m = H ° »<*-*•>• 

a(t) 

= a 0 exp (<f>l ^ 2 (t)) = a 0 exp (^ 0 2 # 2 (t))] , 

(3) 


respectively, where M p i is the Planck mass and to is arbitrary epoch when a = a 0 , <f> = <f > o and 
H = H 0 . The above expression is a good approximation when <j) satisfies the following inequality. 

r)*^)**-*- < 4 > 

For <f> > fa the evolution of (f> is dominated by quantum fluctuations [18] and at <f> = <f>^ the time 
variation rate of <p or $/ $ becomes as large as the expansion rate H so that inflationary expansion 
terminates. The present horizon scale left the Hubble radius when (f> « 3 M v i corresponding to 
the e-folding number of inflation after this epoch to be about n « 60. In order for the density 
fluctuations to be small enough on this scale m should satisfy the constraint m < 10 13 GeV [19]. 

As is seen in (2), for H(t) > m /^/ 3 or <p > M p ij\f^K time variation rate of the Hubble 
parameter is so small that the evolution of the scale factor is indistinguishable from that in de 
Sitter spacetime during each span of several expansion times. Hence we may utilize various results 
of quantum field theory in de Sitter spacetime. 

Here we consider the evolution of a scalar field x with the Lagrangian 

•c* = j(flx)’ - v bc], V{x) = jtx 1 - »’) J + (5) 

with a positive coupling parameter ( to the scalar curvature 72.. Depending on the number of 
components of 3i this system allows a domain wall (j = 1), a string (j = 2), a monopole 
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(j = 3), or a texture (j = 4) solution if £72. C Au 2 . In discussing the phase transition, however, 
we concentrate on the single- component case. Generalization to multi-components cases will be 
discussed in §6. 

In the inflationary stage the scalar curvature 7 Z is given by 

n = 121T 2 + 6H 2 - 2m 2 . (6) 

The symmetry of x is restored if £ R > An 2 and second-order phase transition takes place as 7 Z 
decreases gradually. 


3 Evolution of the scalar fields 


In this section, we trace the evolution of x to discuss when its sign in each domain becomes 
fixed so that we may predict where topological defects will appear after completion of the phase 
transition. The symmetric state x = 0 becomes classically unstable when £72. becomes smaller 
than An 3 . Evolution of the scalar field after this epoch can be divided to two stages. In the first 
stage when its potential at the origin is still nearly flat, quantum fluctuations govern its evolution 
and its amplitude grows gradually. In the second stage typical amplitude of x becomes so large 
that its dynamics may be determined classically and its fate in each domain becomes predictable. 
We may regard that x' s sign is fixed in most domains at this epoch and follow its evolution 
classically thereafter. 

In terms of the mean square value of a scalar field $ with a constant mass M [20], 


<*’(*)> = 


3 H* 
8t r 3 M 3 

*L t 

4t r 3 ’ 

3 H* 


r, ( 

2 M 2 t\ 

1 — exp 1 

3 H ) 


exp 




for M 2 > 0, 

for M 2 = 0, 
for M 2 < 0, 


( 7 ) 


{ 8tt 2 |M 2 | 

it has been naively concluded in the previous literatures [13, 21] that the first stage lasts while 
|lkf 2 eff | < H 2 , where M 2 eff = V"[x = 0] = — Au 2 + ^72. is the mass squared at the origin, because the 
linear growth of mean square field with time is characteristic of the era dominated by long-wave 
quantum fluctuations. However, this criterion is not valid in the present case in which M 2 eff is 
time-dependent. 
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In order to find a more appropriate way to estimate when the evolution of y becomes deter- 
ministic, let us focus on its long-wavelength part which is responsible for the symmetry breaking 
and satisfies the Langevin equation 


du 


Hx] + /(«) 


( 8 ) 


3 H 2 {u) H{u)' 

under the assumption of slow rolling, which is justified if \M^ ceS \ 9 H 2 [23]. Here we have used 

a new time variable u = lna(£) — lna(<o) which is the e-folding number of cosmic expansion from 
t 0 to t. The Hubble parameter squared is given by H 2 (u) = H 2 { 0) — \m 2 u = H% — | m 2 u. In the 
right-hand-side (RHS) of (8), the first term represents classical potential force, while the second 
term stands for random quantum noise with the correlation (f(ui)f(u 2 )) = — u 2 ). The 

above Langevin equation yields the following equation of motion of (x 2 (ti)) 


d , a/ .. 2 MLg , 2/ H 2 (u) 


(9) 


where the term proportional to (% 4 (u)) has been neglected, which is shown to be a good approx- 
imation at the end of this section. Again the first term of the RHS is classical potential force 
and the second term represents quantum diffusion without which (9) would be the same as the 
classical equation of motion of y 2 . Thus by comparing magnitudes of the two terms, one may 
conclude if the evolution of y is dominated by potential force or fluctuations in typical domains. 

Equation (9) may readily be solved: 


<x 2 (-)> = 




-+2«+l 


du' 


r 8 « u l - 


2 m 2 u 

~3 W 


-2r+ a < 


( 10 ) 


If 2 m 2 u <C 3jHq and m 2 At; 2 , 


( 2 m 2 u\ 

V wJ 


-+a<+i 


exp 


f 2Xv 2 \ 
\ 3Hq / 


( 11 ) 


and (10) reads 

<X 3 («)> 


(X 2 (0)) exp 
1 


2(12 £H 2 - Xv 2 ) 




u 


3H* 


4tt 2 2(\2(,Hq - At; 2 ) 


{- 


exp 


2(12£H 2 - Xv 2 ) 


3H 2 


"]} 


( 12 ) 


5 


t 


Since 12 (Hi — An 2 is equal to the mass squared of x, = 0], the above result is equivalent 

to that with a constant mass M xe g[u = 0] calculated by the one-loop field theoretic method in 
de Sitter spacetime (7) 1 . With the above approximation (11), (x 2 ( u )) may increase or decrease 
depending on the sign of 12£Hq — Xv 2 or the choice of t 0 . 

We should hence consider higher order terms in order to discuss transition to the second stage. 
Let us first consider the last exponent of (10), 

2mM - ^ +2 * 




/ 2 m 2 u\ 

V wj 


= G(u). 


We can expand In G(u) as 


In G(u) = —8 £u + 


Xv 2 + 2 *m 2 


m A 


2m 2 1 

(2m 2 \ 2 1 . 

(2m 2 V 

3H 2U+ 2 

V3ifo / + 3 1 

v 3-^0 / + "; 


2(Av 2 + 2*m 2 — 12£Hl) 2m 2 (Xv 2 + 2 £m 2 ) 2 


3 m 


9ifn 4 


Av 2 + 2*m 2 / / 2m 2 

(iSf 


3m 2 


(13) 


The approximation (11) corresponds to adopting only the first term of (13). However, in the case 
M 2 eS [u = 0] = 12^LTq — 2*m 2 — Xv 2 is small, in which we are interested, the second term may 
dominate the first and G(u) grows. For simplicity, let us choose to so that M 2 efr = 0 or 

(Av 2 + 2*m 2 ) 


H l = 


12 * 


at that time. Then for small u, G(u) is given by 


G{u) S£ exp 


32 ( 2 m 2 


[An 2 + 2*m 2 

Thus G(u ) starts to grow exponentially at u = uj with 


u 


(14) 


(15) 


u f = 


Au 2 


32* 2 m 2 


32* 2 


An 2 


c = 


m 


(16) 


Similarly the integrand of (10), which is denoted F(u'), may be approximated as 


P / ir /N 2m 2 (m 2 + 2*m 2 + Au 2 ) 

' nf( “ ) 


u + 


3ff 2 (Xv 2 + m 2 + 2 (m 2 - 12*ff 2 ) 


12 


2(m 2 + 2 *m 2 + A v 2 )m 2 


1 The equivalence between stochastic inflation method and fleld theoretic method has been discussed in [24]. 
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(Av 2 + m 3 + 2 *m 2 — 12^o ) 2 Av 2 + m 3 + 2 (m 3 ( ( 2 m 3 

— — °IUfff 


2 (m 2 + 2 *to 2 + A v 2 )m 3 

32* 2 to 2 ( 1 \ 2 m 3 

V + 8?j + 2Av 2 ’ 


3m 2 


Av 2 


( 17 ) 


where in the last expression we have used our assumption of An 2 ^ m 3 . Thus the integral yields 


J^F(u')du' « | 


u 

Uf 


u ^ «/> 
u > u h 


(18) 


The next task is to estimate (x 2 (0))> which may depend on the initial state of the universe in 
general. One may expect, however, that the amplitude of (\ 2 ) is 0((H/2tt) 3 ) at the epoch when 
M 2 eff = if 2 , independent of the initial condition, if the inflationary expansion lasts long enough 
by that time. We have numerically calculated (x J (0)) using (10) starting with this condition for 
various values of c = A v 3 /m 3 ranging from 5 to 200 and * = 3/16. (The reason why we adopted 
3/16 as * will become evident in the next section). As a result the inequality 

< x » (0) > < S.M 

tx w> ~ 4lr ,y m2 . 

was always satisfied. Hence (x 2 (0)) is Unimportant in (10) for u > Uf and we may conclude 


(x’M) « 5 


An 2 


exp 


32* 2 7n 2 
, Xv 3 


-u 


(19) 


47T 2 y 32 * 2 m 2 

Inserting this into (9), we find that classical potential force begins to dominate quantum fluctuation 
at tx fa Uf. 

The condition u = uj is realized at 

'c + 2( 


t — to + 


771 


12* 




when the Hubble parameter, 




6? 


is related with M 3 ^ as 


|M|1| - UV2( Sc ( ) 


(20) 
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One can estimate the e-folding number of inflation after t = tf as 



c_ _ 1 __ J_ fc 
4 4£V2' 


( 21 ) 


In order that topological defects, which are thus produced, may leave observable traces, n/ e should 
be smaller than » 60, otherwise they are inflated away from the present horizon. Hence this sets 
an upper bound on c. On the other hand, it should be at least larger than zero for the above 
arguments to be valid. Thus c should satisfy the following inequality for £ = 3/16. 


3 < c < 90. 


( 22 ) 


So far we have entirely neglected contribution of { x 2 { u )) to M 2 ff and the term proportional 
to (x A ( u )) = 0(( X \ U )) 2 ) in (9)- It is justified if 


w 2/ w A 2 u 2 /Aw 2 
(X («)) ~ 1 28 v /27r 2 ^ 3 V m 2 


< m 2 , 


(23) 


which we may regard as a constraint on A: 


A < mViir’fVi. 


(24) 


It is easy to satisfy the above inequality. Even if we take a rather large value of c, say c = 90, it 
only requires that A < 10" 2 when £ = 3/16. 


4 Power spectrum of fluctuations 

Had there been no correlation beyond the Hubble horizon, we could have identified the correlation 
length of the phase transition with the Hubble length at t = tf. However, the scalar field is 
correlated on various scales due to the inflationary expansion, which is characterized by the power 
spectrum |x*| 2 of long-wave quantum fluctuations given in terms of mode function x*(t). Here we 
calculate it. 

We decompose the operator x as 

X(x, t) = j ^^[akX*(0 e,kX + 4x*(0 e_ikX ]> ( 25 ) 
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where a k and a\ are annihilation and creation operators, respectively, and the mode function 
satisfies the following equation of motion. 

Xk{t) + 3 Hxu{t) + X*(0 + M xeffXfc(0 = 0. (26) 

Using a new variable Xk = &^(t)xk the above equation is rewritten as 


d 2 


+ 


(m)’ 


3 • 




Xk = o. 


(27) 


In the inflationary stage, when slow rolling condition |if| H 2 is applicable, we may regard 
the Hubble parameter to be constant in several expansion time scales during which physical 
wave number k/a{t) decreases exponentially. Moreover since we are interested in the regime 
|M 2 eff | < H 2 where quantum fluctuations are important, equation (27) can be approximately 
solved by the Hankel function of rank 3/2 and the positive frequency mode corresponding to that 
in the Minkowski vacuum for k Ha(t) reads 

**<•>“ Gsbj) ’ 

where the mode function Xk(t) is appropriately normalized according to 


( 28 ) 


»( 0 *I (0 - xlWxkW = *■ 


(29) 


While the above expression (28) is a good approximation, the comoving mode k can shift from 
short-wavelength regime to long- wavelength regime. In the latter phase, in which (A:/a(t)) 2 is 

negligible in (27), one may solve it by means of the WKB method [22]. 

1 ^ 1 
Xk{t ) = A k ai (t k ) ex P [J t S(t')dt' 4 - BuaUtiA ( 


+ B k at(t k ) ex P ~J t S{t')dt' , (30) 


where 


S(t) 




+ 




(31) 


9 H 2 3H 2 

The WKB approximation is valid provided |5| S 2 . In the present model 5(t) is explicitly 

written as 
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In the case £ = 3/16 = 0.1875, which is slightly larger than the conformally coupled case £ = 
1/6 = 0.166..., S(t) is time-independent: 


S = J Xv 2 — -m 2 . 


(33) 


This is why we have taken ( = 3/16 in §3 and, for simplicity, we will mainly study this case in 
the following as well. 

The coefficients Ak and Bk may be determined by the initial condition at t = t k and normal- 
ization condition. One may determine the coefficient of the leading term Ak with the help of (28), 
which is explicitly written as 








- s 3 S/ . 


(34) 


where the last expression is the approximate solution for the case the physical wavelength of k- 
mode is larger than the horizon scale ( k < aH ). With the approximation 5(t) — §-BT(t), (30) 
reads 

Xk(t) = Ak<£ (t fc ) exp H(t - t fc )j = Aka* (t). (35) 

Thus from (34) and (35) we can determine Ak as 


A h = 


EM 

2 k 3 ’ 


(36) 


where the unimportant phase factor has been omitted. Thus the mode function in the long- 
wavelength regime is given as 




(37) 


One may get this result independently of the choice of the connecting time £* as far as both (28) 
and (30) are reasonable approximations at t = £*. Hereafter we take £* to be the epoch fc-mode 
leaves the de Sitter horizon, that is, 


k = H{t k )a{tk) = H k a k . 


(38) 
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For the case £ = 3/16, (37) reads 

*(‘) = { -4^ [*2 - *’«] + T [ff ‘ * * (i)l ) ' (39) 

This mode function begins growing at 

= H c (40) 

at which epoch M^ eff = - n j-. Because of the gradual decrease of the Hubble parameter, it is 
slightly after x = 0 becomes unstable that x* starts to grow. 

The amplitude of (39) in the case £ = 3/16 and c = 5 is shown in Fig. 1, which depicts relative 
amplitude to that of a massless minimally-coupled scalar field in the exact de Sitter spacetime 
given by (Hf/2k 3 )^. 

5 Distribution of the topological defects 

5.1 Method and results of numerical simulations 

Having estimated when evolution of x becomes dominated by the potential force in §3 and calcu- 
lated the power spectrum of its fluctuations at an arbitrary time in §4, we are now in a position 
to apply them to find the spectrum of topological defects produced. Following our two step ap- 
proximation described in §3, this is accomplished by calculating spatial distribution of the scalar 
field at t = tf. Thanks to the properties of the vacuum state realized as a result of inflationary 
expansion, x( x >t/) is classically given by summing up its Fourier modes as 

XCM/) = H^XkCM/) + C.C. ; £xk(x,*/) = |Sxk(*/)!e ikx+ ‘ Vk , (41) 

k 

where |£xk(t/)| is a random value which has a Gaussian distribution with dispersion |x*(t/)| 2 
is a random phase-factor. 

If we find the sign of x positive (negative) at a certain point at t = tf, it will fall down to the 
minimum of x = + u (— v) as phase transition proceeds, since its dynamics has been deterministic 
by this time in typical domains. Thus we may consider that domain walls are produced between 
two regions with opposite signs of x at this time. 
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We have first done three-dimensional simulations as illustrated in figure 2A, in which we have 
calculated the value of x(x,t/) through (41) at each point of 32 3 lattices using the power spectrum 
of Xk with £ = 3/16 and c = 5. We adopted fast Fourier- transform (FFT) method to speed up the 
calculation [25] and summed up 32 3 modes of 8x k with k ranging from 27T.H//32 to 27riT/, where 
H J X is the separation of neighboring lattice-points taken to be the Hubble length at t = tj. In 
the figure dotted points will fall in the plus minimum or x = v and empty points in the minus 
minimum after the phase transition so that domain walls will be present between them. As is seen 
there, there are structures on various scales. This is in contrast to figure 2B, in which the sign of 
X has been randomly assigned to each lattice-point corresponding to the case of a thermal phase 
transition with the correlation length Hj 1 . 

In both cases, however, simulation boxes are dominated by walls with infinitely large surface 
area. This is simply because both states of x = are realized with equal probability and hence 
it is unlikely that walls axe so distributed that islands of plus minimum exist in the sea of minus 
minimum or vice versa. Thus it is not appropriate to examine the area-distribution of walls in 
order to discuss the character of phase transitions. 

Hence we instead focus on distribution of separation of each wall. For this purpose we do not 
have to calculate the value of (41) at all points in three dimensional space. But all we should find 
is its value along a line, which enables us to employ a much larger simulation box. In fact we have 
used a box with (2 13 ) 3 lattice points and examined the distribution of domains along a line with 
2 13 points. We have done simulations for the following two cases: 

• A. Power spectrum is given by (39) with £ = 3/16 and c = 5. 

• B. A scale-invariant power spectrum, 

2 = s 

2k 3 ’ 

is assigned corresponding to the case with a massless scalar field in de Sitter spacetime. 

In both cases we have summed 2 13 x 100 = 819200 independent mode functions in (41) with k 
ranging from 2 ttJT// 2 13 to 2i rHf. Simulations have also been done in the case a box with (2 12 ) 3 
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points. As a result no artificial boundary effect was observed. In addition simulations in which 
the upper bound of k is taken to be irHf have been carried out and no effect of small-scale cutoff 
emerged. For comparison we have also studied the case 

• C. Sign of x is assigned at random at each point. 

We have done simulations for 500 times for each case and their results are depicted in Figure 3(A- 
C). We can obtain various consequences from the graph. As is seen there domain-size distribution 
in phase transition during inflation (A and B) is again very different from that in thermal phase 
transition (C). In the latter case we can fit the distribution with an exponential function 

n(s,tf)ds oc exp(— 0.67s)ds, (42) 

where n(a,tf) is the number density of domains with separation a ~ a + da. This is simply because 
the number of domains with size aHj 1 is proportional to the probability of having the same sign 
of x for a times in succession: 

n(s,tf)da oc da oc exp(— ln2 • a)da oc exp(— 0.67a)da. 

Thus our numerical simulation agrees with the analytic estimate. 

On the other hand, in the case B we may fit the distribution with a power-law 

n(a,tf)da oc a~ p da, with p = 3.0. (43) 

We can extend the above one-dimensional result into three dimension as follows. First let us 
define the effective correlation volume of a domain between two walls with separation a by V = a 3 . 
Then, since the probability to find a structure with this correlation volume by our simulation is 
proportional to its effective surface area or a 2 , number density of domains with correlation volume 
V ~ V + dV, n(V,tf), is related with n(s,t/) by 

n(a,t f )da oc 

3 

Thus we have 

n(V,t f )dV oc n(s,tf)s 2 ds oc n(s,t/)d 3 s oc — — . (44) 

s 
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Therefore in the case scale-invariant fluctuation is assigned, we find that the spectrum of the 
resultant topological defects is also scale-invariant. Analytic arguments to derive this distribution 
will be given in the next subsection. 

Finally the plot A in Fig. 3, for which more realistic spectrum has been assigned, has almost 
the same profile as the plot B for a < 80 but larger domains are suppressed because amplitudes 
of fluctuations with larger wavelength are smaller. From Figs. 1 and 3 we may conclude the 
fluctuation mode with \J k 3 \xk{t f)\* smaller than ~ 10% of its maximum is unimportant and the 
spectrum of topological defects on such a large scale should approach that predicted by white- noise 
fluctuations. 


5.2 Analytic interpretation of a scale-invariant spectrum 


Here let us consider how the above scale-invariant feature of topological defects is explained 
analytically. To make the following discussion simple and clear let us focus on a patch of space 
which was homogeneous with x = 0 at t — 0. Let us further assume that the space may be 
divided into many domains of horizon volume H~ 3 with different values of x a t each time. Then 
probability distribution function of x is given by 


P [xOM) = x] = / . e x P 

J2nG(Q,t) 


X 

2G(0,t) 


G(0,t) = (x 2 (x,t)) = — Ft. (45) 


Furthermore joint probability distribution reads 


= Xi;x(x;,<) = Xj] = 


1 ex f Xi +x 2 - M r ^)xo 

),t)J 1 — p 2 (r,t ) 1 2G(0,$)[1 — p 2 (r, t) 


2irG(0,t) v 'l — p 2 (r,t) l 2U^,t)[l - p\r,t)] ) 

where p(r,t) = G(r,t)/G(0,t ) = (x(xi,t)x(xj,t))/G(0,t) with |x* - Xjj = r. The condition of 
existence of a wall at a certain point is that x changes its sign there, whose probability, F(t), is 


given by 


F(t ) = / d Xi dx;2P[x(xi,t) = Xi;x(*j>0 = Xjb |x< - x,| = H l , 

J — oo J0 

1 [G(fl'- 1 ,t)1 _ 1 / IN 1 /T 

7T G(0,t) ir \ HtJ 7 r V HF 


= i arccos % a. I arccos (l - ±) Z U* (47) 

7T G(0,t) 7T V HtJ 7T V Ht V ’ 

where the last approximation is justified if Ht 1. Thus F[t) depends on t only weakly even 

though background space is expanding exponentially. 
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By assumption of homogeneity at t = 0, the largest possible wall at time t has a physical size 
of a = H~ 1 e Bt , so that F(t) is related with ra(s,t) as 

r£T- 1 e H * 


F(t) w j n(s,t)a 3 da, 


(48) 


where t — dependence of n(s,t) reflects the fact that quantum fluctuations are generated continu- 
ously to produce and destroy topological defects successively in the course of cosmic expansion. 
Since this is a rather stationary process without any exponential instability, n(a,t ) should not 
depend on t exponentially just as F(t) does not. Thus from (48) we expect that n(a,t ) should 
be proportional to a -3 to eliminate the exponential dependence on t of the upper bound of the 
integral. In this way we conclude that distribution of walls is scale-invariant and given by 


for Ht > 1. 


«(s, t)a 2 da « 


1 d?a 
(. Ht Jtl 3 "’ 


(49) 


6 Discussion and Conclusions 

We have studied generic features of a second-order phase transition during inflation which is a 
more natural and attractive scenario to produce topological defects at GUT scale than the Kibble 
mechanism. We have considered a specific model in which \ i s coupled with spacetime curvature 
72. and phase transition takes place due to gradual decrease of the Hubble parameter. In this model 
we have clarified the nature of long-wavelength quantum fluctuations by explicitly calculating the 
power spectrum and (x 2 )* 

In terms of the latter quantity we determined when classical potential force surpasses quantum 
fluctuations in typical domains and investigated the distribution of the field at that epoch through 
numerical calculations, by which we examined distribution of topological defects. 

As a result we have found that they have an almost scale-invariant distribution with small- 
and large-scale cutoffs. The former is naturally given by the comoving scale corresponding to the 
Hubble radius at the epoch time evolution of x becomes dominated by classical potential force, 
while our simulations have shown that the latter corresponds to the scale on which the amplitude 
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V 0b 3 | x *(*/)| 2 of fluctuations becomes smaller than ~ 10% of its maximum value as is seen in Figs. 
1 and 3. On the larger scale the fluctuation amplitude becomes so small that the distribution 
should approach white-noise spectrum. 

Note that though Hodges and Primack reported similar results recently, i.e., scale-invariant 
distribution of topological defects produced during inflation, their treatment is not satisfactory in 
two aspects. First they have assumed that topological defects with a certain comoving scale are 
produced only when that scale leaves the Hubble radius during inflation. Contrary to their claim, 
however, defects are produced and destroyed on various scales continuously through successive 
generation of fluctuations as discussed in §5. Second they have not given cutoffs of the distribution 
explicitly. As is seen above we have significantly improved these points through our analyses. 

Finally we stress that our basic results are also applicable to other topological defects such as 
strings or textures in which we are more interested. For example, strings are loci on which both 
real and imaginary parts of string-forming complex scalar field changes sign. Hence if self coupling 
A of the field is small enough we can identify strings with intersections of two different kinds of 
independent “domain walls.” Thus with the same reasoning as §5.2 they also have a scale-invariant 
distribution with cutoffs. Furthermore, since the model parameters must satisfy 

Xv 3 


m 


3 ~ 


< 90 


(50) 


in order that strings thus produced are not diluted too much by subsequent inflation, we have 


A S 10 ' 4 (iraHlS^v)’ 


(51) 


so that the condition of A being small is automatically satisfied for cosmologically important 
strings. 

Thus unlike the conventional scenario, strings or textures in the present model have almost 
scale-invariant spectrum from the beginning without any evolution but can serve as .seeds for 
galaxy formation as well under natural choices of model parameters. 
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Figure Captions 


Fig. 1 Magnitude of the mode function |x*(t/)i for £ = ^ and c = 5 as expressed in equation 
(39). Plotted values are relative amplitude to (.By /2A: 3 )?. Abscissa corresponds to the wavelength 
2 tt a(tf)/k in unit of Hj l . 

Fig. 2A Distribution of domains in plus vacuum. Dotted points depict regions with x = + v - 
Values of the field are determined by giving correlations specified by (39). The box size is 32 3 and 
lattice separation is equal to Sj l . 

Fig. 2B Same as Figure 2A but the sign of x each point is given randomly. 

Fig. 3 Distribution of the separation of domain walls. This figure shows the sum of 500 times 
simulations of the length 2 13 . Plotted values are integrated one in each logarithmic interval. 
A(0) : x(x,£/) is calculated using the power spectrum given by (39). 

B( X ): x( x ,ty) is calculated using the scale-invariant power spectrum. Solid line is the power-law 
fitting whose power index is —2.0 due to integration. 

C(A): The sign of x( x >*/) is assigned randomly with equal probability corresponding to the 
thermal phase transition with correlation length Bj 1 . The plot fits an exponential function 
depicted by the dashed line. 
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